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We present a method to calculate the soft function in Soft-Colhnear Effective Theory to NLO 
for A'^-jet events, defined with respect to arbitrarily complicated observables and algorithms, using 
a subtraction-based method. We show that at one loop the singularity structure of all observ- 
able/algorithm combinations can be classified as one of two types. Type 1 jets include jets defined 
with inclusive algorithms for which a jet shape is measured. Type II jets include jets found with 
exclusive algorithms, as well as jets for which only the direction and energy are measured. Cross 
sections that are inclusive over a certain region of phase space, such as the forward region at a 
hadron collider, are examples of Type II jets. We show that for a large class of measurements the 
required subtractions are already known analytically, including traditional jet shape measurements 
at hadron colliders. We demonstrate our method by calculating the soft functions for the case of 
jets defined in rj-cj) space with an out-of-jet pr cut and a rapidity cut on the jets, as well as for the 
case of 1-jettiness. 



1. INTRODUCTION 

It is well known that the perturbative expansion of jet 
cross sections generally contains large logarithmic terms. 
At each order in perturbation theory, there are powers 
of logarithms of ratios of scales, such as the ratio of the 
jet mass over the jet energy or the ratio of the invariant 
mass between two jets over the energies of the jets. In 
many cases, the presence of these logarithms spoils the 
convergence of perturbation theory, and for that reason 
such large logarithms are often resummed to all orders in 
the perturbative expansion. 

Several techniques have been developed to allow 
the resummation of these logarithms [T], and recently 
it has been shown how soft-coUinear effective theory 
(SCET) [2H3 can be used to resum logs. The first step 
in resummation is to factorize an A'^-jet cross section into 
hard, jet, and soft functions 

dCTAT - BnHnX [J ]" ® S'jv ■ (1) 

Here Bn denotes the Born-level cross section in full QCD, 
Hn reproduces the virtual corrections of full QCD, while 
the J's and Sn together encode real emission diagrams 
in the coUinear and soft limits. There is one jet function 
for each jet in the final state, and both the jet and soft 
functions depend on the algorithm used to define the jets, 
as well as the observables measured. Note that Eq. ([T]) 
must be modified in the case of hadron collision, since 
there will also be PDFs and, for some measurements, 
beam functions [6j. At tree level, the hard, jet, and soft 
functions are trivial, but each has to be calculated order 
by order in perturbation theory. 

The simplest jet definition involves exactly two jets, 
each consisting of all particles in one of the two hemi- 
spheres defined by a plane perpendicular to the thrust 
axis, and is typically only used in e"'"e~ collisions. In this 
case resummation has been achieved at NNNLL |X, . 

For more complicated jet definitions, however, the re- 
quired calculations are more involved, and in many cases 



we do not know the NLO results for the jet and soft func- 
tions. One counter example is cone or inclusive fc^-type 
algorithms in e+e^ collisions, where the distance mea- 
sure is the angle 9 of each particle with respect to the 
jet axis, and the total energy outside of the jets is less 
than A. For this example, the jet and soft functions are 
known to ©(a^) [8]. 

Unlike for e+e", jet definitions at hadron colliders are 
usually required to be boost invariant. As a result, dis- 
tance is usually measured in r]-(j) space, where rj and 
are defined with respect to the beam axis, and there is a 
restriction on the total |pt| outside of jets, rather than 
energy. Other, more inclusive jet definitions, such as A'^- 
jettiness [S], also have complicated dependence on the 
kinematics of the event, such that the NLO results for 
the jet and soft functions are not known. The absence 
of these results has been one of the biggest hurdles in 
deriving more precise predictions for jet cross sections at 
hadron colliders, and the solution to this problem is the 
topic of the current paper. 

A generic A^-jet cross section is defined by a jet algo- 
rithm, which identifies the regions of phase space belong- 
ing to each jet and includes restrictions on the out-of-jet 
radiation, and possibly one or more jet shapes, which 
measure functions of the final state particles in each of 
these jets. For an A^-jet cross section, where an observ- 
able is measured for m jets, the 0{as) contribution to 
the soft function can be written as 

^T.I^M^^M;k), (2) 

where A^^- is defined as 

A^,, = -5V''T,.T, 27r,5(fc^)g(fc°) . (3) 
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In Eq. ([2]), the sum goes over all pairs of soft emission 
sources i and j, which includes the jets and possibly the 
beams, in the case of hadron collisions. The function 
encodes the action of the jet algorithm A, represents 
the measurement M, and ct^ is the value of the z*'' jet 
shape. We will restrict ourselves to measurements of the 
form 
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where restricts the soft gluon k to be part of jet i, 
forces the k to be outside of all jets, while still con- 
tributing to the jet cross section (typically enforcing a 
maximum energy or px value). The function mea- 
sures an observable in jet i, while setting the observables 
in the other jets to zero 

A^{k)^S{ak~ak{k))l[S{ai) . (5) 

Note that for a jet where no jet shape is measured, = 

As an example of a jet algorithm, consider a cone algo- 
rithm of size R at an e'^e~ collider, where the out-of-jet 
energy is restricted to be less than A. This gives 
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l-EQco„e(fc)W(fc°<A). (6) 
^ k=l / 



At one loop, it can be shown that the soft function is 
at most divergent, which means that the result can 
be written as 



5« (^,X;K„})-(aO)n^K) 



(7) 
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where we have defined the moments of the soft function 
with respect to the crfe as 

, ^0 



d'^fc 
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Here we have assumed that all ai are normalized to 1. 
Note that we are suppressing the dependence of the mo- 
ments on the algorithm and measurement. In terms of 



the integral in Eq. ([2]) , the zeroth moment can be written 
as 

d-^fc 
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where the subscript "meas" denotes the set of m jets for 
which an observable is measured. The higher moments 
for measured jets (fc G meas) are 
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E 



d'^k 



(2vr) 



^%(fc)a^'(fc)e^(A:)0K(A:)<l), 



(10) 



while for k ^ meas, (cr^>°) = 0. 

Note that the moments (ct") are in general divergent, 
with the divergences arising from the soft (/c° — > 0), 
collinear {ui ■ k 0), and ultraviolet (fc° 00) lim- 
its. This means that they have to be calculated an- 
alytically, with regulators for all divergences, which is 
in general only feasible for the simplest jet algorithms. 
However, consider the difference of two soft functions, 
S{A,M; {<7m}) and S{A,M;{am}), which are defined 
with different algorithms and measurements. We will 
look at the difference 



-Dk = i'Tk] 



(11) 



where (erf) are moments of 5^^^ and (<t") are moments of 
S''^\ Note that in practice it is often useful to define the 
moments cr" and differences for each color structure 
separately, which we will denote by [(t"]^ and 

In the soft limit, fc° — 0, IR safety dictates that 
CTi(fc) — > 0. This implies that all higher moments (and 
therefore their differences) vanish. For the zeroth mo- 
ment, the observable theta functions are trivially satis- 
fied, such that 
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= 0. (12) 



Here we have used the fact that gluon will not be vetoed 
as its energy goes to 0. This implies that it is assigned 
either to be in a jet or to the out-of-jet region, such that 

J2iLo ®A = 1 lo"" algorithm. 

In the collinear limit, n^-fc — )> 0, IR safety once again 
tells us that a'i{k) — 0. In addition, for fixed fc°, any 
IR-safe algorithm will assign the emission to jet i in this 
limit, such that ~ 5ik- This leads to the desired 
result 



lim 

^i■k^O 
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(13) 
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The final limit, — )■ oo, requires a more careful anal- 
ysis. As the energy of the gluon goes to infinity, any out- 
of-jet cut will veto the emission, which means we only 
need to consider radiation in the jets. The difference of 
the zcroth moments is 



Jhn^(/0(fc)-/0(fc))cx [e.A(fc)-©l(fc)] (14) 

fc^mcas 

+ E [<dAik)0{<Jk{k) <l)~Q^{k)e{ak{k) <l) 

A; tineas 

while for the higher moments we find 

hm (x<f(fc)ej^(fc)0K(A:) < 1) (15) 

fc'J— >-oo 

-5I'(fc)e_|(fc)0(5,(fc)<i). 

In contrast to the soft and collinear limits, where ob- 
servable dependence is trivial, the UV divergences are 
naively sensitive to both the algorithm and observable 
definitions. However, we find that there are two possible 
cases for jets: Type I, where the bound on the observable 
is more restrictive 

Typel : lim Q\{k)e{(Ji{k) < 1) = lim e{ai{k) < 1) , 

(16) 

and Type II, where the jet algorithm is more restrictive 

Typell: lim e^{k)e{(7i{k) < I) = lim e_^(A:) . 

(17) 

Here we have defined the limit of the algorithm restric- 
tion and the observable in the UV. For jets where an 
observable is not measured, the jet algorithm is, by def- 
inition, more restrictive, such that these jets are always 
Type II. 

Clearly, if S and S do not agree on which jets are Type 
I versus Type II, the difference of moments will not go 
to and 2? will be divergent as fc" — > c». Assuming they 
agree, Eq. ([T4|) can be simplified, giving 



lim {l"{k)~P{k)) cx 
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(18) 



pSTypo II 

while Eq. (jlsj) gives 
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eTypc I) 



(19) 



lim 



aUk)0{<Jk{k) < 1) -a^^{k)9iak{k) < 1) 



Here we have used the fact that, for Type II jets, all 
higher moments are [TUj. We now see that, as fco ~^ co, 
if limfcQ^oo o-fc = limfco^oo for ah Type I jets and 



limfeo^oo 6^ = hmfeo^oo ©_| for all Type II jets, the in- 
tegrand vanishes as k'^ ^ oo. 

These results can be used to calculate a desired soft 
function numerically, given an analytically known soft 
function. Combining Eqs. ([t]), ([s]), and (11) gives 
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where we have used that the (log(Ti/(Ti)+ distribution 
must agree [lOj . The notation am Cm means that, in 
Ui are to be replaced by Ui for Type I jets, while 
for Type H jets, Y{^<i,^K^i) becomes Ilj^mcas '^(^O- 
As long as the subtraction soft function has the same 
UV observable dependence as the target for all Type I 
jets, as well as the same UV algorithm dependence for 
all Type II jets, the I?" are finite and can be calculated 
numerically. If S is known analytically, the result for S 
can be fully computed at one loop. We will now illus- 
trate this procedure with several examples. For a longer 
discussion, see |10) . 

For our subtraction function, we will use the results 
from [8]. By combining the pieces of the soft function 
calculated therein, we can construct a subtraction for 
any soft function at NLO that meets the following two 
conditions: observables for Type I jets must be symmet- 
ric about the jet axis and behave like an angularity |llj 
with a < 1 in the UV, and all Type II jets must be found 
using either a cone or inclusive fcy-type algorithm, with 
the 0-angle relative the jet axis as a measure. 
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FIG. 1: The coefficient [D%j normahzed to - for 

three equally-spaced jets as a function of the subtraction jet 
size Rs. 

Using the notation of [S], our subtraction function is 
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Here, m is set to be the number of Type I jets in S , while 
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the are chosen such that \imka^rx,{o'k —Ta) = for all 
Type I jets. 

As an initial check of our method, we seek to repro- 
duce the known analytical result for of three equally 
separated jets, each with energy 125 GeV, found using 
a cone algorithm of size R measured as the angle with 
respect to the jet and an out-of-jet energy cut of A = 50 
GeV. Jet thrust (a = 0) will be measured on each jet. 
Our subtraction will use size Rs, the same out-of-jet en- 
ergy cut, and /i = A = 50 GeV. As we will show in more 
detail in [TO], the result reduces to a simple numerical 
integral over the solid angle. We see in Fig. [T] that the 
known result is reproduced, within numerical error due 
to Monte Carlo integration. 

In Fig. [2j we show the result of of the soft function for 
a single measured jet with a = at a hadron collider. 
The jet is defined with an rj-cf) algorithm, and has energy 
125 GeV. For simplicity we will choose the beams to have 
energy of 125 GeV as well. To veto other jets, we use a 
Pt cut of 20 GeV for \r]\ < 5, while remaining inclusive 
for |?7| > 5. We show three separate curves for R = 
0.4, 0.7, 1.0. The plot is shown as a function of the angle 
of the jet with respect to the beams. Note that to remain 
inclusive for \r]\ > 5 requires the use of Type II jets in 
the subtraction. 
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FIG. 2: The soft function for an rj-cp jet algorithm as a func- 
tion of the angle of the measured jet with different R for the 
case when the emitters (ij) are the two beam directions. 

Using the same kinematic configuration as in the rj-cj) 
example (two fixed back-to-back directions, with a third 
direction of varying angle) , we show in Fig. the result 



for 1-jettiness [S]. Since there is a measurement in each 
region, only Type I jets are used in the subtraction. 

We have seen that, at one loop, subtractions can be 
constructed that allow for new soft functions to be cal- 
culated numerically, using previously derived analytical 
results. Using only calculations that exist in the litera- 
ture, the soft function can be calculated for a number of 
phenomenologically relevant observables and algorithms, 
including jet shapes found using r]-(/) algorithms and A'^- 
jettiness, for arbitrary N. In order to calculate the soft 
function for most current measurements at hardon collid- 
ers, both Type I and Type II jets are required. We should 
also note that a similar argument can be used to calculate 
jet functions in the presence of algorithms; however, this 
is generally less relevant as algorithm corrections tend to 
be power suppressed. 

12 h 
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- (23) 



0.2 0.4 0.6 O.S 1.0 12 1.4 jet 

FIG. 3: The coefficient V° for one-jettiness, shown as a func- 
tion of the angle of the measured jet. The labels refer to the 
three (ij) emitters, with 1 and 2 referring to the beams and 
3 to the jet. 
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